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I. INTRODUCTION 

In the past two decades much progress has been made 
in our understanding of the QCD phase diagram. This 
has been achieved in a combination of first principle con- 
tinuum computations in QCD, see e.g. PUS], lattice sim- 
ulations, see e.g. [6-8J, as well as computations in low en- 
ergy effective models for QCD, see e.g. [9HT4]. All these 
methods have their specific strengths, and the respective 
investigations are complementary. In some applications, 
the methods have even been applied in combination. 

In recent years the first principle continuum approach 
as well as the model computations have become far more 
quantitative. This is the more interesting as at present 
they are one of the methods of choice for finite density 
and non-equilibrium investigations that are necessary for 
an understanding of heavy ion collisions. A full quanti- 
tative understanding of the phase diagram of QCD from 
these methods requires a good quantitative grip on the 
thermodynamics and finite temperature dependence of 
Yang-Mills theory, for a review on thermal gluons see |15j . 
for applications of Dyson- Schwinger equations (DSEs) for 
Yang-Mills propagators see e.g. [T6| fTT]. 
However, quantitatively reliable continuum results for 
correlation functions in Yang-Mills theory at general tem- 
peratures are presently not available. This concerns 
in particular the transition region of the confinement- 
deconfinement phase transition, occurring at the critical 
temperature T c . The same applies to thermodynamic 
quantities such as the pressure. Computations in Yang- 
Mills theory done with hard-thermal loop or 2PI resum- 
mations do not match the lattice results for temperatures 
T < 3T C , see e.g. [18] and references therein. In turn, 
they are reliable for T ^> T c . 

In the present work we put forward an approach with 
the functional renormalisation group (FRG), for QCD- 
related reviews see [TJ [T9 25j . The FRG incorporates 
non-perturbative effects by a successive integration of 
fluctuations related to a given momentum scale, also 
at non- vanishing temperature, for applications in Yang- 
Mills theory see e.g. pU I26H29] . It is applicable to all 



temperatures, in particular also for temperatures T < 
3T C . The present implementation of the FRG incorpo- 
rates only thermal fluctuations, see e.g. [29J, and utilises 
the vacuum physics at vanishing temperature as an in- 
put. This also allows us to study the effect of global gauge 
fixing issues in the Landau gauge due to the Gribov prob- 
lem, for detailed discussions see [3Q| I5T]. A summary of 
the present work and results is also presented in |32j . 

In Section [IT] we give a brief introduction to the FRG- 
approach to Yang-Mills theory. In particular, it is illus- 
trated at the example of the two-point correlation func- 
tions, how flow equations for correlation functions are 
derived from the flow of the effective action. 



Section [TTT] is devoted to the important issue of locality 
of an FRG-flow. We show how to minimise the momen- 
tum transfer present in the diagram, and hence to min- 
imise the sensitivity to the approximation at hand. We 
also discuss the factorisation of diagrams at large exter- 
nal momenta that is relevant for subleading momentum 
corrections. 

In Section IIVI we introduce the flow for thermal fluc- 
tuations as that of the difference of flows at finite tem- 
perature and vanishing temperature, see [29]. The local- 
ity introduced in Section [TTT] is used for improving the 
stability of this difference, as well as having access to a 
trivial initial condition for the thermal flow at large cut- 
off scales. Variants of which have been very successfully 
used already for the physics of scalar theories as well as 
fermion-boson mixtures, for reviews see [24, 33, 34j. 

The approximation used in the present work is detailed 
in Section [V| and some computational details are dis- 
cussed in Section [VI} Results for propagators and vertices 



are presented in the Sections 
to lattice results 



shortly summarised in Section 
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VII| They are compared 
Our observations are 



II. YANG-MILLS CORRELATION FUNCTIONS 
FROM RG-FLOWS 

In this section we introduce the functional RG setting 
for Yang-Mills theories which we use for the computation 







FIG. 2: Functional flow for the effective action. Lines with 
filled circles denote fully dressed field dependent propagators 
([7]). Crossed circles denote the regulator insertion d t Rk- 



FIG. 1: Regulator R(p 2 ) and its (logarithmic) cut-off scale 
derivative d t R(p 2 ) = R(p 2 ). 



of propagators and vertices. Central to this approach 
is the scale-dependent effective action T^, where quan- 
tum fluctuations with momenta p 2 > k 2 are already in- 
tegrated out, and k denotes an infrared cut-off scale. By 
varying the cut-off scale k one thereby interpolates be- 
tween the classical action in the ultraviolet and the full 
quantum effective action in the infrared, where the cut- 
off is removed. An infinitesimal change of k is described 
by a flow equation for T^, and the interpolation can be 
achieved by integrating this flow. 

Within the standard Faddeev-Popov gauge fixing pro- 
cedure in the Landau gauge, 



0, 



1 N 2 - 1 

5 ' * ' 5 C ' 



(1) 



the classical gauge fixed action for SU(N C ) Yang-Mills 
theory in Euclidean space-time is given by 



S- 



4/ f!U ± ILV 



c a d^D ab c b 



(2) 



with f x = J d 4 x. In ^ we have already introduced the 
ghost action with the ghost fields c a and c a . The field 
strength F and the covariant derivative D are given by 



pa 



— r) A a 



acb 



9f 



abc a b 



A A c 



= S ab d^gf acb A 



(3) 



The gauge fields are transveral due to the Landau gauge 
([T]). Hence operators and the inverse of operators are 
defined on the transversal subspace. The suppression of 
the infrared fluctuations is achieved via a modification of 
the propagator with S — » S + ASk and 



AS k 



a a -nab \b 



> b c b . 



(4) 



with f 



J d A p/(2it) A . The momentum-dependent regu- 
implement the infrared cut- 
« k 2 for gluon and ghost re- 
spectively. They carry a Lorentz and gauge group tensor 
structure and are proportional to a dimensionless shape 
function r, schematically we have R(p 2 ) oc p 2 r(p 2 /k 2 ). 



lator functions R ab v and R ab 
off at the momentum scale p 2 



An example for R(jp 2 ) and its logarithmic scale deriva- 
tive is given in Fig. [I] As we work in Landau gauge, ([!]), 
we can restrict ourselves to a transversal gluon regulator, 

v»RfM = 0. 

With the notation ip = (A, c, c) we write the flow equa- 
tion for the Yang-Mills effective action r&[A, c, c] at finite 
temperature T as 

Wkiv] = \i GfM{v,p)d t R%{p) 
V 

-£ G ab [^](p 1 p)d t R ba (p) 1 (5) 

where t = ln/c, and the momentum integration measure 
at finite temperature is given by 



d 3 p 



with po = lnTn , 



(6) 



where the integration over po turns into a sum over Mat- 
subara frequencies. Both, gluons and ghosts have pe- 
riodic boundary conditions, ip(xo + 1/T, x) — (p(xo,x), 
which is reflected in the Matsubara modes 2nTn with a 
thermal zero mode for n = 0. At vanishing temperature 
we have J - ^ I - The dependent propagator for 

a propagation from tpi to y>2 is given by 



1 



M- (7) 



In Q we also used the regulator function in field space, 
Rk }( p 1( p 2 with 



r 



k,A a A h u 



r>ab 



R 



k,c a c b 



-R 



k,c b c a 



R ab 



(8) 



The flow is finite in both, the infrared and the ultravi- 
olet, by construction. Effectively, the momentum inte- 
gration in ([5| only receives contributions for momenta 
in the vicinity of p 2 < k 2 . Consequently, it has a re- 
markable numerical stability. The flow solely depends on 
dressed vertices and propagators, leading to consistent 
RG-scaling on either side of (5). Diagrammatically, the 
flow in ([5| is depicted in Fig. 2 

The structure of the flow ([5|, or Fig.|2j is more appar- 



or 
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FIG. 3: Field derivative of the propagator. The solid line 
stands for either a gluon or a ghost propagator and the vertex 
is depicted by a filled circle. 



ent in a formulation with the fields cp. Written in these 
fields the cut-off term reads ASk[^p] = \y> • Rk ' where 
the dot indicates the contraction over species of fields in- 
cluding space-time or momenta, Lorentz and gauge group 
indices. Note that the cut-off term now includes two 
times the ghost term which cancels the global factor 1/2 
and leads to Q. In this short hand notation the flow 
equation takes the simple and concise form 



d t T k [<p] = -TrG[<p](p,p)-d t R k (p), 



(9) 



where the trace includes a relative minus sign for the 
ghosts. Eq. ([9| is well-suited for structural considera- 
tions. 

Flow equations for propagators and vertices are ob- 
tained from ([5]) or ^ by taking derivatives with respect 
to cp. The computation of ip- derivatives on the rhs of ([9| 
only requires the equation for the field derivative of the 
propagator and the notion of the full one-particle irre- 
ducible (1PI) n-point functions. The latter are defined 

by 



,Pn) 



6<p(pi) ■ ■ ■ <p(p n ) 



The field derivative of the propagator is given by 



Sip 



• G[<p] 



(10) 



(11) 



and diagrammatically depicted in Fig. [3j With the above 
equations one readily derives the flow of 1PI n-point func- 
tions. Again, this can be nicely illustrated diagrammat- 
ically at the example of the propagators, which are the 
key objects in the approach. We take two field derivatives 
of the flow equation. Applying this to its diagrammatical 
form displayed in Fig. [2] we are led to the flow equations 
given in Fig. [4j The flow equations only contain one loop 
diagrams in the full propagators, which is a consequence 
of the exact one-loop form of the flow for the effective 
action, see ([5| and Fig. [2] Moreover, they only depend 
on the full vertices ( flQ| ) . Analogously to the flow of the 
two-point functions, the flows of n-point functions only 
contain one-loop diagrams that depend on full propaga- 
tors and vertices. 

The propagators and vertices are purely transver- 
sal due to the transversality of the gauge field. The 



FIG. 4: Flow equations for the ghost and gluon propagators. 
Vertices with filled circles denote fully dressed vertices. All 
propagators are fully dressed, the filled circles for the internal 
ones have been omitted for the sake of clarity. Crossed circles 
denote the regulator insertion dtRk- 



purely transversal correlation functions vanish by con- 
tracting one of the Lorentz indices with its momentum, 

te)^ r m- T /iz---/i™ = for z = l,...,m. Note that in 
general n ^ m due to the ghosts. Since the propa- 
gators are purely transversal in the Landau gauge, the 
purely transversal correlation functions in the Landau 
gauge form a closed system of flow equations: the flow 
of a purely transversal correlation function only depends 
on purely transversal correlation functions and carry the 
whole dynamics. Any observables can be built-up from 
the purely transversal correlation functions. 

In turn, the flow of correlation functions with at 
least one longitudinal direction, r^ n ^ L , depends on both, 
r( n ) T and r( n ) L . Moreover, they also obey modified 
Slavnov- Taylor identities, see [2T| [22j [45j 06] and refer- 
ences therein. This situation is summarised in 



dtT (n)T 

<9 t r (n) L 



(p)„r 



(n)L 



Fiow^[{r^}], 
Fiow^[{r^ T ,r^ L }], 
mSTi n [{r^ T ,r^ L }], (12) 



where m < n + 2. The modified Slavnov- Taylor identities 
converge to the standard ones for vanishing cut-off scale 
k = 0. The above system comprises the full information 
about the correlation functions in the Landau gauge. In- 
terestingly, the hierarchy of flow equations for the purely 
transversal correlation functions can be solved indepen- 
dently and carries the full dynamics of Yang-Mills theory. 
In this context we remark that vertex constructions aided 
by Slavnov- Taylor identities implicitly utilise an assumed 
uniformity of correlation functions, that is 



d p r^ < 00 



for all (pi 



,Pn) 



(13) 



This works well in perturbation theory but has to be 
taken with a grain of salt in the non-perturbative regime. 



It is also worth emphasising that (12) does not depend 
on the way how the Landau gauge is introduced. The 
Landau gauge can be also represented as the limit of 
covariant gauges with the gauge action l/(2£) j^d^A^) 2 
and gauge fixing parameter £. Here £ = signals the 
Landau gauge. In this case in general one also introduces 
a regularisation for the gauge mode, schematically given 



3 



by 



R L 



1 



lim -U^p z r(p z /k z ), 



(14) 



with the projection operators 



\1V 

n^(p) 



PuPv/p 2 ■ 



(15) 



on transversal and longitudinal degrees of freedom, re- 
spectively. Still, the longitudinal mode does not play any 
role for the flow of correlation functions as G • dtR gauge ■ 
G —> for £ — )> and lim^ G is purely transversal. 
Note however, that l/2TiR gauge ■ G does not vanish. 
Upon t-integration it gives the thermal pressure for the 
gauge mode, namely the Stefan-Boltzmann pressure of 
N 2 - 1 fields, see [25]. 



III. LOCALITY OF RG-FLOWS 

The loops on the rhs of Fig. [4] only received contri- 
butions from momentum fluctuations with p 2 < k 2 due 
to the regulator insertion dtR, see Fig. [I] However, the 
external momenta are not limited by such a constraint. 
Indeed, for large external momenta, p 2 /k 2 ^> 1 the flow 
factorises at leading order: the tadpole diagram with the 
four-gluon vertex tends to a constant. The other tadpole 
diagrams vanish with powers < k 2 /p 2 . The related four- 
point functions are not present on the classical level and 
hence decay at least with k 2 /p 2 . This intuitive statement 
can be proven easily with the help of the respective flows 
and the factorisation present there. Here, we explain the 
factorisation at the relevant example of the three-point 
function diagrams, see also Fig. [5] for the example of the 
flow of the ghost propagator. These diagrams factorise 
for large external momenta, one factor being the uncut- 
ted internal line evaluated at p 2 . For p 2 /k 2 — >• oo we 
have, 



Tr 



(GRG) (q) • r< 3 > (q, p + q)-G(p + q)- r< 3 > (q + p, q) 

(grg) (g)-r( 3 )(o,p)-G(p)-r( 3 )( P ,o)' 



Tr 



+ higher order terms , 



(16) 



where the trace also integrates over loop momenta q. In 
(16) we have assumed that all momentum components 



are suppressed by the regulator: R is a function of p 
and not e.g. of solely spatial momentum squared, p 2 . 
Note also that for kinetic and/or symmetry reasons the 
leading order in line two of (16) may vanish. This even 



supports the factorisation. The interchange of integra- 
tion and limit and hence the factorisation works as the 
diagram is still finite with the uncutted line being re- 
moved. 



The general factorisation in (16) can be nicely illus- 




p/k > 1 




X 



p 



FIG. 5: Factorisation in leading order for large momenta for 
the first diagram (cutted gluon line) in the flow of the (inverse) 
ghost propagator in Fig. [4] . The triangle stands for the prod- 
uct of the two vertices at q — and reads —p^p^'f acd f bc d . 



the ghost propagator Fig. [4] In Fig. [5]we display the 
factorisation of the first diagram in FigT^J The terms in 
the second line of (16) or on the rhs of Fig. [5] are already 
subleading as the flow usually is peaked at momentum 
scales p 2 < k 2 . However, they are already quantitatively 
relevant at vanishing temperature but turn out to be cru- 
cial for the correct thermodynamics, in particular for the 
slow approach to the Stefan-Boltzmann limit for large 
temperatures. Potentially, they also play an important 
role for the thermodynamics in non-relativistic systems, 
where they supposedly relate to the Tan-relations [47U49] 
in the context of many-body physics; for FRG-reviews see 
e.g. [31134]. 

Still we have to reconcile the above polynomial de- 
cay with external momenta p with the well-known ex- 
ponential decay of thermal fluctuations with the stan- 
dard suppression factor exp(— m/T) in the presence of a 
mass scale m. In the present case, this mass scale can 
be either the cut-off scale k or the physical mass scale of 
Yang-Mills theory, Aqcd, which is directly linked to the 
critical temperature T c of the confinement-deconfinement 
phase transition. The exponential thermal damping 
factor originates from the full Matsubara sum, and is 
strictly not present if only the lowest Matsubara fre- 
quencies are taken into account. A four-dimensional reg- 
ulator depending on four-dimensional loop momentum 
q 2 = (27rTn) 2 + q 2 cuts the Matsubara sum, and hence 



trated diagrammatically at the example of the flow of 



only leads to the polynomial decay of the flow. The expo- 
nential suppression is then built-up successively with the 
flow. The above properties can be already very clearly 
seen and understood at the example of perturbative one- 
loop flows. 

The full Matsubara sum is reintroduced for regulators 
only depending on p 2 . They are frequently used in finite 
temperature applications of the FRG as they allow for 
an analytic summation of the Matsubara sums if only 
the trivial frequency dependence is taken into account, 
see e.g. [T9} 129} [501 151]. However, in the present work we 
consider full propagators and vertices, that is with non- 
trivial frequency and momentum dependence. Moreover, 
as lower-dimensional regulators introduce an additional 
momentum or frequency transfers in the flow we therefore 
refrain from using them here. 

Apparently, the large momentum contributions also 
weaken the locality of the flow present in the loop mo- 
menta as they induce a momentum transfer: the flows at 



4 



a given cut-off scale k carry physics information about 
larger momentum scales. In turn this entails that any 
local approximation does not fully cover this momentum 
transfer. Now we utilise the freedom of reparameterising 
the effective action as well as the choice of the regula- 
tor in order to reinstall the locality of the flow for the 

(2) 

two-point functions 1^ . This minimises the systematic 
error of a given approximation, see [21J. To that end, we 
rewrite the regulator term in Q as follows, 



r k [(p], see [H], 



2j A (p))R b «(p) 



G ab [4>)(p,p) (d t + 2^ c (p))R ba (p) , (23) 



AS k 



^k,a fiu ^k,b 



C K aR ab C k , h . (17) 



The fields <\> = (Ak, c k ) in ([17]) relate to the cut-off 
independent fields ip = (A, c, c) in the classical action via 
a cut-off and momentum dependent rescaling, 



where G[(j)){p,q) = (lY j [0] + R) (p,q) denotes the full 
regularised propagator for the propagation of (/>, see ^ 



Hp) 



^(pMp) 



with dt<t>(p) = %(p)<t>(.P), (18) 



(2) ~ (2) 

with ip — »• (j) and 1^ —> 1^ . The functional flow in 
(23) looks rather complicated, but it is simply a repa- 
rameterisation of the standard flow in ([5|, or (|9|. In the 
condensed notation introduced in Section [IT] this is more 
apparent, the flow equation (23) then reads 



where the derivative is taken at fixed (p. We will use the 
natural definition % — %- Eq. (18) implies 



d t + <t> • % ■ 



■ T.k I 



= -Tr G[4>] ■ (dt 



2%) ■ R k 



%(p) = dtlogZ^ip) , and R = ZR. 



(19) 



(24) 

Eq. ( 24 ) illustrates that we have only reparametrised the 



The field reparametrisation in (18) does not change the 



effective action, in particular the regulator term ASk 
does not change. It simply amounts to rewriting the ef- 
fective action in terms of the new fields, 



fields in a scale-dependent way. 

Taking two derivatives with respect to <\>\ (p) and </>2 (q) 
of ( 23 ) at vanishing ghost fields and constant gauge field 



(20) 



we schematically get the flow 

(^+^ 1 ( P ))fg^(p) 



F1 owi 2 l(p), 



Then, ^-derivatives f ( k n) of the effective action T k = t k where the rhs of {25} stands for the Mp) and 



are given by ( 10 ) with (p — » (j>, 



(25) 
Hi) 



fi n) (Pl,...,Pn) = 



derivative of the rhs of (|23|, 

5 2 



5<p(pi) ■ ■ ■ 5<j){p n ) ' 



(21) 



As the fields <j) and cp only differ by a momentum depen- 
dent rescaling with Z 1 / 2 , the correlation functions are 
related by a simple rescaling with powers of Z 1 / 2 , 

n 

r^( Pl ,...,p n ) = Y[z^ 2 (p i )r^(p 1 ,..., Pn ). (22) 



The rescaling with Z is up to our disposal, and we shall 
use it to minimise the momentum transfer in the flow 
equation of the two-point function by eliminating the 
subleading terms in the flow exemplified in (16). Note 
that this does not remove the related contributions, the 
momentum transfer is still present but does not feed- 
back directly in the flow, see [21J. This is also elucidated 
below. 



Flow 



and 



(2) 



^TrG[0] • (dt + fj^-Rk) , (26) 



V<f>(p) = 2%(p) = d t log Z k {p) , 



(27) 



is the 'anomalous' dimension of the propagator related to 
the rescaling of the fields with Z . In (I25| we have used 
that for vanishing ghosts the two-point functions are di- 
agonal/symplectic in field space. The only non- vanishing 



components are and f ^] 



^(2) 



The two- 



point functions for constant gauge fields are also diagonal 
in momentum space, that is 



rfW = r^(p)(27r) 4 %-g), 

Flow (2) (p,g) = F\ow {2) (p) (2tt) 4 5 (p-q), (28) 



^(2 



The flow equation for the effective action now receives 
further contributions from the /c-dependence of the fields 
(j). We are finally led to the following flow for T^tM = 



and the relation \22\ for the two-point functions reads 
rW(p) = Z(p)fW(p). (29) 

Now we can come back to the question of locality of the 
flow. Our aim is to remove the large momentum tail 
displayed in ( 16 ) from the flow in order to keep its locality 



5 



in momentum space, i.e. to minimise the momentum 
transfer. This is achieved by demanding 



fcff (p) 



p 2 >(\k) 2 



which implies 



V<t> (?) = 



Flow (2) (p) 



rl 2 V) 



0(p 2 - (Xk) 2 ) . 



(30) 



(31) 



Eq. ( 30 ) entails that the momentum transfer is switched 



off for momenta larger than the cut-off scale. The factor 
A controls this scale and can be used for an error estimate. 

After having computed the localised correlation func- 
tions T^ n \ we can derive the correlation functions r^ n ^ 
via rescaling with powers of Z, see (21 ) ,(|29j). The scaling 



factor is computed by integrating f)k defined in (27) 



{i: 



Z k (p;T) = Z k=0 (p;T) exp\ / dt' fj k > (p; T) > .(32) 



This determines Z k up to a fc-independent function, and 
we use 



Z<t,,k=o(p;T = 0) = I 



(33) 



For the choice (33) the two sets of correlation functions 



agree in the vacuum at k = 0. 

As the flow of general correlation functions can be writ- 
ten down solely in terms of T^ n \ the relation (30) with 
(31) eliminates the momentum transfer (in T^j from 
the flow. Note however, that a remnant of it is still 
present via the factor 2^ = fj^ on the rhs of the flow 
(24). For regulators that decay sufficiently fast for mo- 



menta p 2 ^> k 2 this is quantitatively negligible. Indeed, 
for regulators which vanish identically for momenta big- 
ger than Xk the momentum transfer now is described 
solely by fj^ip) and decouples completely. We also re- 
mark in this context, that the above construction and 
the definition ( 33 ) leading to p0|) with A = 1 can be de- 



duced by evoking functional optimisation for momentum- 
dependent approximations, see [5T]. The above heuristic 
arguments entail that optimisation restores the locality 
of the flow also in general momentum- and frequency- 
dependent approximations. 



IV. THERMAL FLUCTUATIONS 

The flow equation ([5| includes both, quantum as well 
as thermal fluctuations. For the present purpose we 
are interested in the thermal fluctuations, for reviews on 
thermal FRG see e.g. pj2[52j[53]. Thermal fluctuations 
are encoded in the difference of the flows at finite and at 



vanishing temperature, see e.g. 
d t AT T ^}= ±Tr GM'dtR 



-Tr GM'dtR 



where 



Ar T>fc [0] = r fc [0]| T - r fc [0]| 



T=0 



T=0 

(34) 



(35) 



accounts for the difference between the effective action 
at vanishing and at finite temperature. Due to the ther- 
mal exponential suppression the flow (34) should have 
locality properties with respect to the scale k = T. As 
discussed in the previous section, locality is important 
for the quantitative reliability of a given approximation. 
In the present work we implement this idea as follows: 
We use the vacuum physics at vanishing temperature as 
input for the flow. A given set of correlation functions 
rL=o &t T = and k = can be integrated with the flow 
(pi) at vanishing temperature in a given approximation 
up to a large momentum scale k = A, 



r£= (pi,...,Pn) 



flow 



T=0 



r^ A (pi,...,p n ) 



T=0 



(36) 



In the approximation at hand, (36) defines the initial con- 



ditions r 



(n) 



T=0 



which give the correct vacuum corre- 
lation functions if integrated to k = 0. The ultraviolet 
scale A is chosen such that all thermal fluctuations are 
suppressed given the maximal temperature to be consid- 
ered, T max . This implies 



A 



< 1. 



(37) 



Then, switching on T < T max does not change the initial 
conditions at leading order, i.e. 



1 



AT^ A (p u ...,p n )=0 + O 



(38) 



where d n is the canonical dimension of T^ n \ and all mo- 
menta are of order A or bigger, p? > A 2 , see e.g. [54] 
for scalar theories. This suggests that we can flow ATt,/c 
from the trivial initial condition ( 38 ) to vanishing cut-off 
k = 0. Note also that within such an approach to ther- 
mal fluctuations, it is only the difference dt AT^/c which is 
sensitive to the approximation at hand. We illustrate the 
procedure outlined above in a heuristic plot, see Fig. [6] 



In this context it is also worth discussing the rapid- 
ity of AT^ — >• for large cut-off scales k = A. This 
is linked to the questkm of locality raised in the previ- 
ous Section 



III 



In Fig. p|it is indicated that AT^ start 



I 1 i k 

to significantly deviate irom zero at the thermal scale 

k ~ 2ttT. Above and in the previous section we have ar- 
gued that we only have a polynomial suppression for the 



flow AT^. Note also, that the longer AT^ survives at 



(n) 



T,k' 



large scales k — >• A the more sensitive it will be to the 
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< 1 




FIG. 6: Flow at vanishing temperature from k = to k = A, 
and flow at T ^ from /c = A to k = with A/T » 1. 
The flow is described in theory space, and the axes label (or- 
thogonal) couplings/observables Oi which serve as expansion 
coefficients of the effective action, e.g. G\ — r^ 2 \p = 0). The 
flows start to deviate at k w T, see the discussion below. 



the UV scale we have exponentially suppressed thermal 
fluctuations for the two-point function. As we expect 
a polynomial suppression due to the arguments in Sec- 
tion |TIT[ the use of the computational simple initial con- 
dition amounts to a temperature-dependent renormalisa- 
tion of r^ n) (T) at non- vanishing cut-off scale k. For the 



same reason, (40) also leads to Z^oiP^T) ^ 1, and we 



have to compute it from the flow of fjk- We have 



Z k (p;T) = Z A (p;0)exp<^ / dt' i) k (p;T) 



Jo dt> Vk(p;0) e f£ dt' (fj k (p;T)-fj k (p;0)) ^ ( 4 -q 



where we have used (32) in the second line of (41). The 



relation (41) entails that the rescaling factor Z k= o(p;T) 
contains the thermal part of the momentum transfer. 



approximation at hand. Eventually, the polynomial con- 
tributions to the flow integrate up to the standard ther- 
mal exponential suppression, that relates to the thermal 
distribution functions as discussed in Section Hill 

We conclude that we have a polynomial decay in the 
flow, which indeed plays a role for computing quantita- 
tively reliable thermodynamic quantities. Additionally, 
we can utilise the ideas about locality to preserve as much 
as possible of the exponential decay, that stabilises any 
approximation scheme. We turn to the localised version 



of the flow derived in the previous Section III As the 



flow of T k vanishes identically for momenta larger than 
Xk for all temperatures T we arrive at 



Af^ A (p) = + O( e - A / T ), 



(39) 



where we have neglected the backreaction of the poly- 
nomial decay of the thermal corrections of higher cor- 
relation functions. The two-point function is the most 
important quantity if it comes to the computation of 
thermodynamic observables. It is exactly here where the 
locality-preserving flow pays off in quantitative reliabil- 
ity. 

Finally, we are interested in the correlation functions 
r( n ) at k = 0, which have to be computed from the 
r( n ) via the rescaling with powers of Z k =o. At vanish- 
ing temperature we have used the natural normalisation 
Z k =o — 1, see (32) and (33), and the two sets of correla- 
tion functions agreed at k = 0. At finite temperature we 



initiate the flow at k = A with T/A. Here we have (39) 
for the difference of the localised correlation functions, 
and we use the natural definition 



Z A (p;T) :=Z A (p;T = 0). 



(40) 



(2) 

Eq. (|40|) implies that (|39|) also applies to AT^ ; A : at 



V. APPROXIMATION 



In this section we discuss the approximation scheme 
used to capture the thermal fluctuations and thus al- 
lows for a quantitative computation of the temperature- 
dependence of the propagators. 

In general, the structure of the flow of the effective 
action, see Figs. |2|3[ entails that flows of n-point func- 
tions depend explicitly on (n + 2)-point functions. The 
relevant example for the present work is the flow of the 
(inverse) Yang-Mills propagators displayed in Fig. [4] 

The flows depends on rj^ with n < 4. In other words, 
the functional flow in eq. ([5|, if broken up in the flows 
of n-point functions, constitutes an infinite hierarchy 
of coupled integro-differential equations. For computa- 
tional purposes the system must be closed: the set of po- 
tentially contributing operators must be rendered finite, 
such that the relevant physics is kept in the approxima- 
tion. Moreover, the approximation should be subject to 
self-consistency checks that give access to the systematic 
error. 

First, let us describe the approximation put forward 
here for the flow of the propagators in Fig. |4j before we 
put down the explicit parameterisation/approximation in 

terms of rj^ with n < 4: we keep the full propaga- 
tors and work with self-consistent approximations to the 
vertices which respect the renormalisation group prop- 
erties of the vertices. We also use the flow equation for 
the ghost-gluon vertex, evaluated at the symmetric point 
with p 2 = k 2 , i = 1, 2, 3, for the momenta of ghost, anti- 
ghost and gluon, respectively. We have checked the re- 
liability of the RG-consistent ansatz for the three-gluon 
vertex with its flow at the symmetric point p 2 = k 2 . 
The reliability of the RG-consistent four-gluon vertex is 
checked by its DSE computed in |55| . 



7 



A. Two-point functions and their flows 

Now we proceed with the parameterisation of our ap- 
proximation. We concentrate on the localised two-point 
functions T^ 2 \ that of the cut-off independent fields are 
then obtained with (22), (32) and (41). At vanishing 



temperature the propagators are described by one wave- 
function renormalisation each, ZA,k(p) and Z Cj k(p)- At 
non- vanishing temperature we have to take into account 
chromoelectric and chromomagnetic modes, the respec- 
tive projection operators P L and P T , 



p ^v(P0iP) = C 1 - <W C 1 - <M (<W -P^Pv/P 2 ) > 



(42) 



where II ^ is the four-dimensional transversal projection 
operator, see (15). Thus we parameterise the gluon with 



two wave- function renormalisations Z L j T . The ghost has 
only a scalar structure at vanishing and finite tempera- 
ture. The parameterisation of the gluons and ghost is 
given by 



^(2) 
- A,L 

AST 



(Po,p) 



Zl(po,p)p 2 P L (Po,P) , 
Zt(po,P)p 2 P T (po,p) , 
^cbo,p)p 2 , (43) 



where the identity in colour space is suppressed and the 
Z r s are functions of po and p separately. The parame- 
terisation of I^ 2 ) follows from that of inJ43|, and is 
read off the definition of 4>(ip) in (18) and (|29[), 



t al/t - Zl/t(p) z l/t (p) p 2 > i?> ~ z c ( v ) z c ( p ) p 2 . 

(44) 

The flow equations for the two-point functions have been 
given diagrammatically in Fig. [4j Their right hand sides 
depend on the two-point functions, as well as three- and 
four-point functions. In particular, we have tadpole dia- 
grams which depend on the ghost-ghost and ghost-gluon 
scattering vertices r^L and T^ 2c , respectively. These 
vertices vanish classically, and in a first approximation 
one is tempted to drop the related diagrams. 

However, they can be considered in a rather simple 
way, which we use for the flow of the (inverse) ghost 
propagator: we insert the DSE-relations for r^L and 
r^ 2c in the related diagrams. This provides a DSE- 
resummation of the vertices in a given approximation to 
the flow. After some straightforward but tedious alge- 
braic computations, it can be shown that this turns the 
flow equation for the ghost in Fig. [4] into the total t- 
derivative of the DSE-equation for the ghost, see Fig. [7[ 
This is nothing but the statement that a flow equation for 
a correlation function can be seen as the differential form 
of the corresponding DSE for the ghost two-point func- 
tion in the presence of the regulator term; both describe 
the same correlation function. The ghost-DSE and its 
t-derivative is illustrated in Fig. [7[ The derivative acting 



-R): 



FIG. 7: Flow equation of the ghost propagator from the total 
t-derivative of the renormalised ghost-DSE. The cutted line 
stands for the scale derivative acting on the propagator of the 
corresponding field, see (45). The square denotes the scale 



derivative acting on the full vertex, the vertex without filled 
circle denotes the classical vertex. 



on the dressed propagator gives 

d t G [<j>] = -G [</>] • d t (r< 2 > [<j>] + V) • G [0] . (45) 

The derivative of the bare propagator in the DSE van- 
ishes. The DSE-flow is finite by construction, as it can be 
derived from the manifestly finite ghost flow in Fig. [4] by 
inserting the DSEs for the ghost-ghost and ghost-gluon 
scattering vertices. These DSEs are also manifestly finite 
and require no renormalisation. To see this explicitly in 
Fig. [4] we discuss the single terms. First, we note that 
the total t-derivative of the propagators, d t G, shown in 
( |45| ) act as UV-regularisation of the loops. They decay at 
least with G 2 as dtT^ tends towards a constant for large 
momenta. The total t-derivative dtG decays even more 
rapidly with d t R for large momenta due to (30). This 
reflects again the locality implemented by (30) and its 
practical importance. The last diagram is proportional 
to the flow of the ghost-gluon vertex. It will be discussed 



in detail in the Section |VC| and is displayed schemati- 
cally in Fig. [TO] Here it suffices to say that the vertex 
itself is protected from renormalisation and its flow de- 
cays rapidly for large momenta. 

Both flow equations given in Fig. [4] and Fig. [7] for the 
ghost are exact and are related via the above resumma- 
tion procedure. In the present approximation scheme the 
DSE is more favorable as it only depends on the ghost- 
gluon vertex, which we shall resolve via its flow. In turn, 
the four-point functions do not appear explicitly in the 
total derivative of the DSE, and the related contributions 
are absorbed in the diagrams with cutted propagators. 
Notably, the ghost-ghost scattering vertex r^L disap- 
peared completely from the set of flow equations of ghost 
and gluon propagators, whereas the ghost-gluon scatter- 
ing vertex T^ 2c is still present in the flow of the gluon 
two-point function. 

Note that a similar procedure in the flow of the gluon 
two-point function leads to two-loop diagrams, which are 
commonly dropped in DSE-computations. This is one of 
the reasons why we refrain from resumming the related 
diagrams. 
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B. Vertices and their flows 



For the vertices we first introduce a convenient param- 
eterisation that naturally captures the renormalisation 
group behaviour. To that end, we utilise the wave func- 
tion renormalisation of the gluon and the ghost and write 



V// 
c 



b 



FIG. 8: Ghost-gluon vertex. 



f^(p 1 ,...,p n ) = l[zV 2 (p i )T(p 1 ,...,p n ). (46) 



1=1 



The Z-factors are chosen to be proportional to Z, and 
hence carry the RG-scaling of the vertex as well as the 
momentum dependence of the legs: they carry potential 
kinematical singularities, see e.g. [56J. We choose 



Zl/t(p) 



^L/T 



0) p 2 



^L/T 



(q)q 2 



9=0 



where Zj^a is either Zk l an d ^fe,T? depending on the 



projection P T I L defined in (42) on the respective leg of 



In a slight abuse of notation we have introduced 



j. £ . in cl 0115111 auuoc ui n uiatiuii vvc nave i 

Z(k): the Z-factors in (49) are functions of p^ and p 2 , 
which are evaluated at (48). The freezing scale k s cx 



Aqcd in (49) is chosen such that it is bigger than p pea k- 
We point out that k s only defines the parameterisation 
of the ghost-gluon vertex. 



Z c (p) = Z c {p). 



(47) 



However, the Z L / T (p) are frozen for p Ppe&ki 
where p pea k is the potential turning point of the 

(2) 

inverse propagator 1^ in the infrared defined by 
9 P (p 2 Z L/T (p)) p=Ppeak = 0. The turning point p peak de- 
pends on k and tends towards zero for k —> 00. Without 
this additional constraint Z L j T would turn negative for 
small fc, which reflects positivity violation. At finite tem- 
perature the turning points depend on T and differ for 
Z T and Zl, the turning point of the latter tends towards 
zero for T — » 00 due to the Debye screening. We em- 
phasise that this is done simply for convenience in order 
to avoid the splitting of a positive vertex dressing into 
two negative factors. The Z's take into account the RG- 
scaling of the fields, and reflect the gaps present in the 
gluonic degrees of freedom. 

This leaves us with a renormalisation group invariant 
tensor T. It is regular up to logarithms and carries the 
canonical momentum dimension as well as the tensor and 
colour structure. For the flow of the propagators, Fig. |4j 
T has to be computed for the three-gluon and four-gluon 
vertex, the ghost-gluon vertex, as well as for the four- 
ghost and ghost-gluon scattering vertices. The latter 
two, which are absent on the classical level, are treated 
in terms of exact resummations with the help of DSEs. 
Now we utilise the locality of the flow: it only carries mo- 
menta q 2 < k 2 and is peaked at about p 2 ~ k 2 . Hence, 
we approximate the vertices by evaluating them at the 
symmetric point at p 2 = k 2 and vanishing temporal com- 
ponents, 



(pi)o = and 



^2 7 2 
Pi = k . 



(48) 



menta p with (48) and we set 



Then the Z- factors in ( 47 ) can be evaluated at fixed mo- 



Z k ,A = Z A (k)0(k-k s ) 
Zk,c = Z c (k) , 



■ Z A (k s )6(k s - k) . 



(49) 



C. Ghost-gluon vertex 



It is left to determine the dressings T for the primi- 



tively divergent vertices T^l , T^i 



(3) 

and r^ c , which have 



a classical counterpart. We restrict ourselves to the clas- 
sical vertex structure. The evaluation at the symmetric 
point (48) leaves us with /c-dependent dressing functions. 



For the ghost-gluon vertex we are led to 



T-X^p) 



%k,cAc 



-[Stl{q,p)\t c 



ibe 



z k,cAc ^Qfif Q 

(50) 

where g is the classical coupling, S^ c is the classical 
ghost-gluon vertex derived from (|2|, and q are the 
ghost and anti-ghost momenta respectively, see Fig. [8] 
The /c-dependent factor Zk,cAc is RG-invariant and de- 
fines a running coupling 



a s (k) 



^k,cAc 
47T 



(51) 



with running momentum scale k. If expanded in powers 
of the coupling for large momenta, a s has the one and 
two- loop universal coefficients of the /3-function of Yang- 
Mills theory, where we have used that Z^ y A — > Zk,A for 
large cut-off scales. 
The flow of 

ZqAc is extracted from that of the ghost- 
gluon vertex. Here this flow is computed within a DSE- 
resummation similar to the derivation made for the ghost 
propagator. Again, the DSE-resummed flow is finite 
as it is derived from the standard flow equation for 
the ghost-gluon vertex which is finite by construction. 
Note however, that already the DSE in Fig. [9] is fi- 
nite without renormalisation procedure due to the non- 
renormalisation theorem for the ghost-gluon vertex. Due 
to simple kinematical reasons, it is also present in ap- 
proximation schemes that respect the kinematical sym- 
metries. 

The above arguments allow us to start straightaway 
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FIG. 9: DSE for the ghost-gluon vertex. The vertex without 
filled circle denotes the classical vertex. 




with the simple ghost-gluon DSE, see Fig. [9j which only 
contains one- loop terms. Similarly to the flow for the 
ghost propagator, we turn the DSE in Fig. [9] into a flow 
equation by taking the t-derivative of Fig. [9J The run- 
ning of the vertex allows for a temperature-dependent 
computation of the vertex, which will have an important 
effect in the calculation of the propagators. 

It remains to project the vertex flow onto that of the 
renormalisation group invariant dressing function Zk^Ac 
The projection on the classical vertex structure is done 
by 



FcAc(Ps) 



l 1 cAclfi 



is, 



cAc\jJL 



tq(3) ]de/r Q (3) idef 
A D cAc\v \. D cAc\v , 



p2=q2 = (p-\-q)2= p 2 



(52) 

with an evaluation at the symmetric point at the momen- 



turn scale p s . For the classical vertex r~^ c = S^( c derived 
from (J2|, the dressing is simply unity, Tk,cAc(Ps) = 1. In 
the present approximation we evaluate the vertices at the 
moment scale k, and hence we define 



(3) 



(53) 



Note that the lhs depends on k via the evaluation at 
p s = k but also due to the implicit dependence of the 
vertex on the cut-off scale. The full vertex dressing in 
(53) also includes the dressing of the legs as split off in 
(46). Hence, the dressing function Zk^Ac is given by 



%k,cAc 



1 



- 7 l/2 
J k,A 



ZkA 



k,cAc ' 



(54) 



The flow dfZk^Ac is determined from (54) and is directly 



related to that of the ghost-gluon vertex. Taking the 
t-derivative of (54) leaves us with 



r\ , ld t Z k ,A , d t Z k A QT , 1 

Ot + -— = 1 = Zk,cAc = Ctl fc,cA 



2 Z K 



ZkA 



The scale-derivative of the full dressing V^Ac is propor- 
tional to the flow of the vertex but also to the derivative 
with respect to the momentum at the symmetric point, 

d t T k ,cAc = [dt^cAciPs) +Psdp s FcAc(Ps)] ps=k • 

Upon integration the flow (|54|) gives us the vertex dress- 



FIG. 10: Flow equation for the ghost-gluon vertex from the 
total t-derivative of the ghost-gluon vertex DSE in Fig. [9] 
The cutted line stands for the scale derivative acting on the 
propagator of the corresponding field, see ( [45] ) . The square 
denotes the scale derivative acting on the full vertex, the ver- 
tex without a filled circle denotes the classical vertex. 



ing of the ghost-gluon vertex at a given cut-off scale /c, 



%k,cAc — z k=0,cAc 



r k dk' 

—rrOt'Zk' ,cAc- 



(55) 



Jo 



For the thermal flows in the present work the initial con- 
dition Zk=o,cAc\ T={d is required. It is here where we take 
advantage of the progress made over the last two decades 
in our understanding of Landau gauge QCD at vanishing 
temperature, for a review see [30]. At vanishing temper- 
ature a one-parameter family of solutions with infrared 
enhanced ghost propagators and gapped gluon propaga- 
tors in Landau gauge has been found. All solutions have 
a gluon propagator with a mass gap m g i UO n oc Aqcd- 
They only differ from each other in the deep infrared for 
momenta p 2 <C Aqcd- There the gluon propagator is 
described by 



Z A (p < Aqcd) oc c(p)- 



gluon 
p 2 



(56) 



where c(p) > 1 is a momentum-dependent function which 
is bounded from below. For all solutions but one c(p) it 
is also bounded from above, these solutions are called de- 
coupling solutions, as the gluon decouples in that momen- 
tum regime. There is one distinguished member of this 
family where c(p) diverges with p 2 + 2KA with ka < — 1, 
see [30J. This solution is called scaling solution as the in- 
frared propagators and vertices are uniquely determined 
by scaling laws up to constant prefactors, see e.g. [3Q| f56l- 

EQ. 

Most importantly, the flow of the ghost-gluon vertex 
in Fig. [lO] is not sensitive to the infrared behaviour of the 



gluon propagator. All diagrams in Fig. 10 vanish in the 
limit k — » with powers of k 2 /m 2 ap , 



lim d t z k ,cAc 



k^O 



OC 



rn 



(57) 



gap 



This also entails that the vertex dressing tends toward 
a constant in the infrared. Hence, we infer that the 
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infrared value of the ghost-gluon vertex is the same 
for the whole class of solutions up to subleading terms 
in RG-transfor mat ions. This can be used to compute 
Zk=o, cAc | T=0 for the whole class of solutions. It has 
been shown that the scaling solution for constant ghost- 
gluon dressings is determined in the FRG up to an RG- 
constant, see [62 J, 



6 k=0,cAc\T = Q 



Za(p)Z?(p) 

z a , s (p)zUp) 



Ana 



s,IR •> 



(58a) 



p=0 



with the scaling wave function renormalisations Za/c,s 
for ghost and gluon propagators, respectively. The cou- 
pling a s ,iR is analytically known, see [58J, 



4tt 2T(-2k)T(k 



i)r(« + 3) n c 



3 (r(-^)) 2 r(2^-i) 



2.97, (58b) 



where k « 0.595. 

For momenta p ^> Aqcd the Z A /c,s(p) tend towards 
the decoupling solutions up to RG-scalings. Demanding 
equivalence for large momenta fixes the relative ultravi- 
olet renormalisation condition. In summary, this allows 
us to fix z 2 = ~ Ac at vanishing temperature for a given 
set of scaling or decoupling propagators in terms of the 
UV renormalisation condition. 

In the light of the ongoing debate about the infrared 
behaviour of Landau gauge propagators in the vacuum 
(T = 0) it is important to stress the following: first, the 
flow of the vertex function is not sensitive to the differ- 
ences of the momentum behaviour of the propagators in 
the deep infrared, p <C Aqcd? as it is switched off for 
k — »• below Aqcd- Second, the above argument leading 
to (58) only relies on the technical possibility of finding 



initial ultraviolet conditions for the flow in the given ap- 
proximation which flow into the scaling solution in the 
infrared. This is trivially possible, see [30 , 62 , 63J. Then, 
the analytical values of the scaling solution fixes (58) for 
both, scaling and decoupling solutions. This does not re- 
solve the infrared problem in Landau gauge Yang-Mills 
theory closely related to the picture of confinement, as 
well as to the resolution of the Gribov problem in this 
gauge. 

At non-vanishing temperature the coupling is sup- 
pressed below the temperature scale k ~ T, see e.g. [28] . 
Furthermore, we have to distinguish between transversal 
and longitudinal gluon legs. If all Matsubara frequencies 
vanish, po = qo = 0, the longitudinal vertex vanishes and 
the distinction is only relevant for non- vanishing Matsub- 
ara frequencies. There, however, all (hatted) quantities 
quickly tend towards their T = counter parts. There- 
fore, we approximate the longitudinal vertex dressing at 
finite temperature by the transversal one, 



7 L 

^k,cAc 



7 T 

^k,cAc 



(59) 



In summary, we have set up a relatively simple flow for 
the ghost-gluon vertex which already covers the quantita- 



FIG. 11: Infrared dominating diagrams in the flow equations 
for the three-gluon vertex and the four-gluon vertex. 



tive features of the full vertex flow. We have also checked 
our approximation by computing the full vertex flow on 
the basis of the results in the present approximations. 
This also provides us with an estimate of the systematic 
error in the present approximation. A full analysis of this 
will be published elsewhere. 



D. Gluonic vertices 

It remains to determine the purely gluonic vertices. 
They are described within a parameterisation similar to 
(50). We have schematically 



7^3 



1 



z k,A 3 ~ ° A3 •> 

9 



7 A 1 C?(4) 

Z k,A 4 ~0 D A 4 •> 



(60) 



where the vertex dressings zj^a 3 an d £&,a 4 relate directly 
to the ghost-gluon dressing z 2 - Ac = Ana s (k), (51), for 
large cut-off scales k ^> Aqcd or large momentaaue to 
two loop universality. 

Indeed, this reasoning has been validated in [55] with 
DSE-equations for the four-point coupling. Hence, we 
have z kj A 4 — z \ cAc f° r mos t of the momentum regime 
with a potential deviation in the deep infrared for 
k 2 / Aq CD <C 1. We parameterise accordingly 



z k,A 3 = z 3 z k,cAc , «s,A 4 = ^4 z k,cAc 



(61) 



where the z^i = 3,4, are functions which are expected 
to approach unity for k 2 > Aq CD , that is 



^,fc»A QCD 



1 



(62) 



for i = 3,4. Their infrared behaviour is determined by 
the only diagram in the flow that does not depend on the 
gapped gluon propagator, see Fig.[TT] The effect of these 
diagrams on the vertices is determined by two competing 
effects: 

First, these diagrams are suppressed due to their colour 
structures and the related vertex dressings are suppressed 
relative to that of the ghost-gluon vertex. In the case of 
the four-gluon vertex this combinatorial suppression is of 
order 10 2 — 10 3 , which is nicely seen in the solution in 
[55J. For the three-gluon vertex the combinatorial sup- 
pression factor also turns out to be of order 10 2 — 10 3 , 
subject to the chosen momenta, details will be published 
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elsewhere. This factor has been determined within a two- 
dimensional lattice study as 0.017, [64] . 

Second, these diagrams grow strong in comparison to 
the respective diagram in the flow of the ghost-gluon ver- 
tex, which involves one gapped gluonic propagator. In 
the decoupling case the flow of the gluonic vertices for 
k — » is proportional to k° up to logarithms. This leads 
to an effective suppression of diagrams with gluonic ver- 
tices as such diagrams also involve gapped gluon prop- 
agators, and hence a suppression with k 2 /m^ ap . In the 
scaling case the diagrams with gluonic vertices decouple 
with powers of the scaling, which can be seen best in the 
DSE-hierarchy. 

In summary, the above analysis entails that we can 
safely drop the related diagrams for k < Aqcd, and above 
this scale the vertices have a dressing similar to that of 
the ghost-gluon vertex. In the vertices we have factorised 
the multiplicative dressing functions Z 1 / 2 for the respec- 
tive legs that are at our disposal. In turn, the relative 
thermal suppression factor is safely encoded in the ratio 
of the full ghost-gluon vertex dressings at vanishing and 
finite temperature, 



J- k,cAc\T 
^k,cAc \t— 



(63) 



Since the temperature dependent dressing factor ZcAc has 
already been included in the definition (61 ), we deal with 



the reduced ratio of the wave function renormalisations, 



r- cAc (k,T) = 



Z k,A Z k,C 


T 


Z k,A Z k,C 


T=0 



(64) 



leading to the final approximation of the gluonic vertex 
dressings by 



z i {k,T)=r- cAc {k,T) i 



\k,T) = ^ Zl (k,T). 

^k,A 



(65) 

The results we show are achieved within the second set of 
vertex dressings, where we additionally switch-off the glu- 
onic vertices linearly below the temperature scale, which 
implements the additional thermal suppression in the 
infrared. This betters the numerical convergence. In 



Section VII we test the sensitivity of the results to this 
switch-off. 

At non- vanishing temperature the structure functions 
za™ also have to carry the difference between the coupling 
to longitudinal and transversal gluons, which is relevant 
for the second choice in (65). The ratio Za/Za goes to 
Zl/t/Zl/t- The validity of the present approximation 
is then tested by computing the flow of Z43 on the basis 
of our results. Note also that z^^Ac only has transversal 
parts as we evaluate the flows at vanishing Matsubara 
frequency. 



E. Regulators 

For a numerical treatment of finite temperature flow 
equations, exponential regulators 



r m (x) 



„m—l 



(66) 



generally provide good numerical stability. The param- 
eter m controls the sharpness of the regulator, which 
directly relates to the locality of the flow. As already 
pointed out in section [TTTJ locality is a central issue. The 
importance of sufficiently local flows is also seen on the 
level of the regulator, where it turns out that the non- 
locality induced by the exponential regulator (66) with 



m = 1 spoils the stability of the flow. On the other hand, 
regulators with a steep descent lead to a slower conver- 
gence of the regularised thermal propagators and vertices 
towards the vacuum ones at T = 0. Indeed, for sharp cut- 
off functions or non-analytic ones thermal modifications 
are present for arbitrarily large cut-off scales A. This in- 
validates the use of the T = initial conditions at the 
initial UV scale A. Accordingly, we use m = 2 in the 
computation, whose form is shown in Fig. [T] The full 
regulators for the gluon and the ghost are given by 



R T/La » 

R a \p) 



5 ab P^ L Z k , T/L p*r 2 ( P yk 2 ), 
5 ab Z c (0)P 2 r2(p 2 /k 2 ), 



(67) 



We 



with the projection operators P L I T defined in (42) 
have chosen Z c (0) instead of Z c (k) for numerical con- 
venience. The ghost renormalisation function Z c (p;T) 
tends towards zero for small momenta and finite temper- 
atures. It is a balance between ghost propagator, gluon 
propagator and ghost-gluon vertex, which prohibits Z c 
getting negative. This balance is better resolved numeri- 
cally with the prefactor Z c (0) in the regulator R c for the 
ghost. At vanishing temperature we have Z^^l = Zk,T, 
and the gluon regulator Rk,A = R T + R L is proportional 
to the four-dimensional transversal projection operator 
II defined in (42), see also the discussion below Q. 



VI. COMPUTATIONAL DETAILS 

The discussion of the last section leaves us with a cou- 
pled set of partial integro-differential equation depicted 
The initial condition is set in the vacuum, 
at vanishing cut-off scale, k = 0. The vertices 



in Fig. [12 

T = 0, 

are determined by (|58bj) and the relations in the Sec- 
tions V C|VD| and the input gluon and ghost propagators 
are taken from [30J, from which we choose a decoupling 
solution. 

In order to solve the flow equation up to a UV scale 
A we adopt an iteration procedure, which is very stable 
for the direction k —> A. Herein the iteration (i + 1) only 
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FIG. 12: Yang-Mills flows for propagators and ghost-gluon 
vertex in the approximation discussed in Section [V] The flows 
for the ghost propagator and the ghost-gluon vertex are DSE- 
resummed and the ghost-tadpole in the gluon equation is ne- 
glected. 



depends on the given solution (i), 



<n) 

■ k,i+l 



k=0,i 



k dk' 



(rg,Flow^).(68) 



Incrementing the number of iterations until the solution 
is stable under further iteration, i.e. rj^ +1 = up to 
desired accuracy, gives the solution of the flow equation. 



The starting point of the iteration we take as T K k 



(n) 



,(0) 



A- T=0 k=Q . Then we converge within the iterations rapidly 

to the initial condition T^ =0 (p 2 ). 

In contrast to the direction k —> A in the zero tem- 
perature case, the flow from high to low scales at fi- 
nite temperature involves instabilities. This is due to 
the self-regulating nature of the equations for the wave- 
function renormalisation of the ghost in combination with 
the ghost-gluon vertex dressing. Their structure is such 
that if one of the quantities becomes small it stops the 
flow of the other and in the following also its own flow, 
which is what happens at finite temperature. However, 
as soon as one of these quantities happens to be negative 
in an intermediate iteration step, the iteration becomes 
unstable, i.e. each iteration step brings the iteration solu- 
tion further away from the correct solution. The system 
is highly sensitive to this, as one has to resolve very small 
values numerically. Therefore, we pursue a more direct 
strategy to solve the flow, namely an evolution of the flow 
according to a Runge-Kutta solver 



ki 



: Flow 



(n) 



(69) 



from the starting condition kjy = A ^> T to ko = 0. In 
the evolution the system reacts on the balancing effect 



FIG. 13: Longitudinal gluon propagator Gl at different tem- 
peratures as a function of spatial momentum. 



between the ghost propagator and ghost-gluon vertex im- 
mediately, and the purely numerical problem of possibly 
negative values of Z c or ZqAc in the iteration is avoided. 
Note that the result is stable under iteration again, as it 
is the exact solution of the equation. The sensitivity to 
the balancing is still present, but shows up in the form 
of a small evolution step size of \ki-i — ki\ < lOMeV. 



VII. RESULTS FOR PROPAGATORS AND 
VERTICES 

In this section we present the results for the ghost and 
gluon propagators and the ghost-gluon vertex. The tem- 
perature is given in lattice units. As has been argued in 
section the results for the ghost and gluon propaga- 
tors show the typical thermal scale 2ttT. Below this scale 
we have significant temperature effects on the momen- 
tum dependence. In turn, above this scale the tempera- 
ture fluctuations are suppressed and all propagators tend 
towards their vacuum counterparts at vanishing temper- 
ature. This also holds true for the ghost-gluon vertex. 
This supports the self-consistency and stability of the 
thermal flows as discussed in section HVl 

The most significant effect can be seen for the chro- 
moelectric and chromomagnetic gluon propagators, that 
is the components of the propagator longitudinal and 
transversal to the heat bath. The zero mode of the longi- 
tudinal gluon propagator at various temperatures is given 
in Fig. 13 as a function of spatial momentum. For low 
temperatures T < 150 MeV we see an enhancement of 
the longitudinal propagator. Such an enhancement is 
also seen on the lattice, [T5j [171 ISSl [36l HH HI]- It has 
been emphasised on the basis of the FRG that specif- 
ically the propagator of the electric mode should show 
critical behaviour, if computed in the background fields 
that solve the non-perturbative equations of motion [35] . 
However, the significance of the lattice results so far as 
well as quantitative details are not settled yet. For higher 
temperatures the longitudinal propagator is suppressed 



13 




Transversal Propagator Gt 

T = MeV 
T = 100 MeV 
T = 200 MeV 
T = 300 MeV 
T = 500 MeV 
T = 1 GeV 



4.5 - 

4 - 

3.5 - 

3 - 

2.5 - 

2 - 

1.5 - 

1 

0.5 - 

- 







Ghost Propagator G c 






T - 





MeV 




100 


MeV 


T 


200 


MeV 




300 


MeV 




500 


MeV 


T - 


1 


GeV 



1.5 2 
P [GeV] 



FIG. 14: Transversal gluon propagator Gt at different tem- 
peratures as a function of spatial momentum. 



FIG. 15: Ghost propagator G c at different temperatures as a 
function of spatial momentum. 



relative to the gluon propagator at vanishing tempera- 
ture. This is the expected behaviour caused by the Debye 
screening mass due to the thermal screening of the chro- 
moelectric gluon. For asymptotically high temperatures 
T ^>T C the chromoelectric gluon decouples. The onset of 
this behaviour at about T w 200 MeV is earlier as in the 
respective lattice computations fT5} fT7| l35l l36| [4T| [42] 
where the thermal decoupling takes place for tempera- 
tures larger than the critical temperature. 

In order to capture this behaviour we have to extend 
our present truncation with a self-consistent inclusion of 
the Polyakov loop background as well as a better res- 
olution of the purely gluonic vertices for momenta and 
frequencies below Aqcd- Both extensions are under way, 
and the results will be presented elsewhere. In turn, for 
large temperature and momentum scales above Aqcd the 
above lack of quantitative precision at infrared scales is 
irrelevant. Here we see quantitative agreement with the 



lattice results, see Fig. 22 



The transversal mode is not enhanced for small tem- 
peratures, in clear distinction to the longitudinal mode. 
It is monotonously decreased with temperature, see 
Fig. 14 Moreover, it develops a clear peak at about 500 



MeV. This can be linked to positivity violation, which has 
to be present for the transversal mode: in the high tem- 
perature limit it describes the remaining dynamical glu- 
ons of three-dimensional Yang-Mills theory in the Landau 
gauge. The infrared bending is more pronounced as that 
of respective lattice results, and its strength is subject to 
the lack of quantitative precision at these scales. In turn, 
for larger momenta the transversal propagator agrees well 
with the respective lattice propagator. The ghost only 
shows a small temperature dependence, in contradistinc- 
tion to the gluonic propagators. This is fully compati- 
ble with the lattice results, see [T5l [TTj [35 l 156) [4T | [42], 
The temperature dependence is hardly evident on the 
level of the propagator, see Fig. [l5j but can be resolved 
on the level of the wave-function renormalisation, see 
Fig. [16] The wave-function renormalisation is slightly 
suppressed, which corresponds to a successive enhance- 
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FIG. 16: Zero mode of the wave-function renormalisation of 
the ghost Z c at different temperatures as a function of spatial 
momentum. 



ment of the ghost propagator at finite temperature. The 
enhancement of the ghost propagator is potentially self- 
amplifying, as it feeds back into the flow of the ghost two- 
point function, see Fig. 
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This would cause a pole in 
the ghost propagator at some temperature, if the ghost- 
gluon vertex did not change. However, in this case the 
flow of the latter is dominated by the diagrams with two 
ghost propagators, see Fig. [l2j and the ghost-gluon ver- 
tex flows to zero. 

This non-trivial interplay of the flow of the ghost prop- 
agator with the flow of the ghost-gluon vertex leads to 
a self-stabilising system and prevents a further enhance- 
ment of the ghost. This effect is crucial for the stability 
of the solution of the Yang-Mills system at finite tem- 
perature. Indeed, for a constant vertex no solution could 
be obtained for intermediate temperatures T « T c and 
above. Thus we conclude that any reliable truncation 
must comprise direct thermal effects also in n-point func- 
tions with n < 3. 

The self-stabilising property of the Yang-Mills sys- 
tem explained above is clearly seen in the temperature- 
dependence of the ghost-gluon vertex. The vertex is sup- 
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FIG. 17: Dressing function ZcAc of the ghost-gluon vertex at 
different temperatures as a function of spatial momentum. 
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FIG. 18: Comparison of longitudinal propagator with mini- 
mal and maximal strength of the coupling of purely gluonic 
vertices. 



FIG. 19: Comparison of transversal propagator with mini- 
mal and maximal strength of the coupling of purely gluonic 
vertices. 




T = MeV 

with bound (T = 200 MeV) 
no bound (T = 200 MeV) 



1.5 2 2.5 3 

k [GeV] 



3.5 



FIG. 20: Comparison of ghost-gluon vertex with unbounded 
and bounded gluonic flow. 



pressed successively with the temperature, see Fig. [T7| 
which in turn ensures the relatively mild change of the 
ghost propagator. Especially the sharp drop-off of the 
vertex at small scales k accounts for the smallness of the 
thermal fluctuations to the ghost propagator. 

In Section |VD[ we approximated the temperature de- 
pendence of the gluonic vertices by an ansatz that com- 
bined all necessary qualitative properties of the vertices. 
In order to check the sensitivity of the propagators at low 
temperatures to this choice, we compare the results with 
a strong suppression with a computation where the mini- 
mal suppression of the gluonic vertices defined in eq. ( 65 ) 



is implemented. This gives an approximate error band for 
the propagators with respect to the ansatz of the gluonic 
vertices. We find that the main impact of the truncation 
of the gluonic vertices is onto the longitudinal propagator 
below T r 



given in Fig. 18 



We have already mentioned 
that our truncation has to be improved in order to quan- 
titatively capture the deep infrared, and the temperature 
dependence of the gluonic vertices is an important ingre- 
dient. 

In contrast to the sensitivity of the deep infrared be- 



haviour of the longitudinal gluon, the transverse gluon 
propagator, see Fig. [l9j as well as the ghost propagator 
hardly feel the modified infrared behaviour of the gluonic 
n-point functions. This is in agreement with the well- 
known fact that the infrared sector of Yang-Mills theory 
in Landau gauge has ghost dominance (for both scaling 
and decoupling solutions) for any dimension d — 2, 3, 4. 

As the wave- function renormalisation of the ghost with 
the minimal and maximal cut off of gluonic vertices can 
not be distinguished by eye, we desist from illustrating 
it explicitly. In contrast to this, the ghost-gluon vertex 
reflects the switching off of gluonic vertices, see Fig. [20} 
however, only in a region where the gluonic contributions 
to the flow are subleading. 

In summary, the scheme-dependence of the gluonic 
vertices affects the details of the thermal decoupling of 
the longitudinal propagator, but has no effects on the 
transversal propagator and the ghost. The sensitivity of 
the thermal decoupling of the propagator to the ther- 
mal decoupling of the vertices is easily understood. In 
turn, transversal gluon and ghost are the relevant de- 
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FIG. 21: Transversal gluon propagator in comparison with 
lattice results [351142] . The lattice data has been rescaled such 
that the T = propagators match at intermediate momenta 
p > 1 GeV. 



FIG. 22: Longitudinal gluon propagator in comparison with 
lattice results [67]. The lattice data has been rescaled such 
that the T = propagators match at intermediate momenta 
p > 1 GeV. 



grees of freedom for asymptotically high temperatures 
where we are led to a three-dimensional confining theory. 
This property can be seen nicely on the lattice, where 
only the transversal gluon is sensitive to the removal of 
confining configurations, see e.g. [65]. We conclude that 
our approximations do not affect the confining physics. 

In the following, we compare the propagators above 
with lattice results [35] [42j [66] . For this purpose, we scale 
the lattice data such that the lattice propagators at van- 
ishing temperature match our normalisation at momenta 
p > 1 GeV. Take notice that we did not use the lattice 
propagator as the initial condition, thus the deep infrared 
of the data deviates from our propagator already at zero 
temperature, which persists also in the propagators at 
finite temperature. Apart from that, there is quanti- 
tative agreement with the lattice data with respect to 
the (temperature dependent) momentum region, where 
the thermal effects appear. In Fig. [21] the transversal 
propagators are compared. The critical temperature in 
the lattice data is T c = 277 MeV. Clearly, we match 
the lattice propagator, except for the strong bending of 
the FRG propagator in the infrared region. This differ- 
ence is a direct consequence of the mismatch of the ghost 
and gluon propagators in the deep infrared at vanishing 
temperature, and relates to the different decoupling so- 
lution chosen here. Apart from the deep infrared, the 
momentum and temperature behaviour of the magnetic 
gluon matches that of the lattice propagator. In con- 
trast to this, the electric gluon on the lattice shows a 
qualitatively different behaviour for temperatures below 
and around the phase transition. Although the longitu- 
dinal propagators agree for T = 0.361T C « lOOMeV, it is 
exactly this region where the uncertainty due to the trun- 



cation for the gluonic vertices is large, as shown in Fig. 18 
for T = 200 MeV. Being aware of a potential truncation 
dependency in the deep infrared of the longitudinal prop- 
agator at low temperatures, we note that in the present 
truncation the electric gluon shows the onset of the en- 



hancement found on the lattice. Increasing the temper- 
ature this feature disappears, and we see a qualitatively 
different effect for temperatures below T c . While the 
continuum result shows a strictly monotonic decreasing 
propagator, the counterpart on the lattice is enhanced in 
the confining regime, but reflects the phase transition in 
form of a rapid decrease at T c . Nevertheless, this deflec- 
tion is expected to be missed in the present truncation, 
as the Polyakov loop potential V(Aq) is pivotal for the 
critical behaviour around the phase transition. In a full 
calculation the inverse propagator is proportional to the 
second derivative of the Polyakov loop T K A J L ~ V"(Aq), 
see [1, 68-70J. This dependence introduces an additional 
screening of the vertex-approximations as well as vertex 
corrections, but in the computations presented here it 
was dropped. This upgrade should also give us access to 
the question of the signatures of criticality in the chromo- 
electric propagator discussed in [15]. Interestingly, such 
a term is absent in the magnetic modes, which do not 
depend strongly on the vertex- approximation. It is sug- 
gestive that the inclusion of Ao will further stabilise our 
computation. However, we defer this interesting compu- 
tation of the propagator in the presence of a non-trivial 
background to future work. 



VIII. SUMMARY AND OUTLOOK 

In the present work we have put forward a functional 
approach for the quantitative computation of full, mo- 
mentum dependent correlation functions at finite temper- 
ature in Yang-Mills theory. This is done within a FRG- 
setting which only incorporates thermal fluctuations. 

We have computed temperature dependent Yang-Mills 
propagators and vertices in the Landau gauge for tem- 
peratures T < 3T C . The chromoelectric propagator shows 
the expected Debye-screening for T > T c in quantitative 
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agreement with the lattice results. For small tempera- 
tures it shows qualitatively the enhancement also seen 
on the lattice [SI H3 E3 ESI EHl HD H2] . However, the 
significance of the lattice results so far as well as quan- 
titative details are not settled yet. The chromomagnetic 
propagator shows the expected thermal scaling and tends 
towards the three-dimensional gluon propagator in quan- 
titative agreement with the lattice. The ghost propaga- 
tor only shows a mild enhancement with temperature in 
agreement with the lattice. In contradistinction we see 
a strong thermal infrared suppression of the ghost-gluon 
vertex which increases with temperature. The results of 
this paper are also summarised in [32J. 

At present we improve on the approximations made 
here and compute thermodynamic observables such as 
the pressure and the scale anomaly. The inclusion of the 



Ao-corrections discussed in Section |VII| should give ac- 
cess to the properties of the chromoelectric propagator 
at criticality, see [35]. Furthermore, we utilise the propa- 



gators obtained here in a combined approach with FRG 
and DSE as well as lattice results in an extension of \71\ . 
This work aims at the QCD phase structure for heavy 
quarks, for lattice results see [6j[72]. Moreover, the glue 
and ghost propagators are a key input for quantitative 
computations in full QCD with 2 and 2+1 flavours at 
finite temperature and density in an extension of [69J. 
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